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We consider the color gauge invariance of a factorized description of the Drell-Yan process 

►^ ' cross section. In particular, we focus on the next-to-leading twist contributions for polarized 

Q \ scattering and on the cross section differential in the transverse momentum of the lepton pair in 

^T' ■ the region where the transverse momentum is small compared to the hard scale. The hadron 

tensor is expressed in terms of manifestly color gauge invariant, nonlocal operator matrix 

Cn , elements and a color gauge invariant treatment of soft gluon poles is given. Also, we clarify 

Cn ■ the discrepancy between two published results for a single spin asymmetry in the Drell-Yan 

cross section. This asymmetry arises if such a soft gluon pole is present in a specific twist-three 

hadronic matrix element. 
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I. INTRODUCTION 



In this article we will address a number of issues concerning the use of nonlocal matrix elements to describe 
Q^ the nonperturbative 'soft' parts that enter in the calculation of hard scattering cross sections. Specifically we 
consider color gauge invariance of these matrix elements in case effects of intrinsic transverse momentum of 
partons are included, necessarily involving nonlocality away from the lightcone. We will study in detail the 
next-to-leading order in an expansion in inverse powers of the hard scale Q for polarized scattering, i.e., order 
l/Q, also referred to as twist-three. This is for instance relevant in the study of the single spin asymmetry 
investigated in Refs. Q and [g|. 

In a factorized description of hadron-hadron scattering processes like the Drell-Yan (DY) process, one can 
5^ ■ work with matrix elements containing nonlocal operators (the local operator product expansion does not apply). 
Cu I At the leading order in inverse powers of the hard scale, operators containing arbitrary numbers of gluons with 
polarization coUincar to the momentum of the parent hadron need to be combined in order to render such 
matrix elements with nonlocal operators color gauge invariant. The gluon contributions sum up to produce 
path-ordered exponentials, also called link operators, with paths running along finite distances in between the 
fields that are situated along a lightlike direction. 

Application of Ward identities to gluons with polarization along their momenta, so-called longitudinally 
polarized gluons, is an important ingredient to arrive at these path-ordered exponentials. At the leading twist, 
but also beyond leading order in inverse powers of the hard scale, one can often perform a so-called coUinear 
expansion of the hard scattering part, such that the partons that need to be considered only have momenta 
collinear to the parent hadron. In that case longitudinally polarized gluons are polarized coUinear to the hadron 
momentum, allowing for straightforward application of Ward identities. However, in some cases (e.g. azimuthal 
asymmetries) one has to deal with partonic transverse momenta and the resulting difference between collinear 
and longitudinal polarization of gluons starts to complicate matters. This difference involves matrix elements 
containing transverse gluon fields and therefore, is only relevant beyond leading order in inverse powers of the 
hard scale. Such a case is for instance the DY cross section differential in the transverse momentum of the lepton 
pair in the region where the transverse momentum is small compared to the hard scale. Despite the mentioned 
complication, we will show for this case that appropriate path-ordered exponentials will result nevertheless, i.e. 
also at the next-to-leading twist (0(1/Q)). 



In case the hadronic matrix elements have a dependence on the transverse momenta of the partons, the 
nonlocaUty of the operators is forced to be off the hght-cone. Consequently, the path-ordered exponentials that 
are needed to render these nonlocal operator matrix elements color gauge invariant, involve paths that are also 
off the light-cone. However, to the order considered here, one finds (as will be shown below) several lightlike 
paths which are transversely separated from each other, depending on the transverse position of the colored 
fields they connect to. The paths extend to lightcone infinity. Such paths are not uncommon, e.g. |§-|§1, but it 
has not been demonstrated before that they also arise in the case(s) under consideration here. 

Hadronic matrix elements involving fields with their arguments at lightcone infinity are usually assumed to 
vanish. In the case of so-called gluonic or soft gluon pole contributions, which are related to transverse gluon 
fields at infinity Ig], this issue needs careful reconsideration. We find that a color gauge invariant treatment of 
soft gluon poles in twist-three hadronic matrix elements 0-0] can be given also. 

In addition to these general results on color gauge invariant matrix elements of nonlocal operators in the 
description of the Drell-Yan cross section differential in the transverse momentum of the lepton pair, we will 
clarify the apparent discrepancy between the results for the power-suppressed single transverse spin azimuthal 
asymmetry in the DY cross section as discussed in Refs. Q and [||. In Ref. [|l| this specific single spin asymmetry 
is assumed to arise from such a soft gluon pole contribution. In Ref. |0] it is shown that the same asymmetry 
can also arise from a so-called time-reversal odd distribution function. It was shown that although time reversal 
symmetry prohibits the appearance of such functions (if one considers the incoming hadron states to be plane 
wave states) the soft gluon pole matrix element can be treated as an effective time-reversal odd distribution 
function, without violating time reversal symmetry. Nevertheless, there remained a discrepancy between the 
results in Refs. |I| and [0 arising from an assumption on the factorized expression for the hadron tensor. We 
will make this assumption explicit and carefully study this issue. First we will discuss the discrepancy between 
the two single spin asymmetry expressions in a bit more detail. 

II. SINGLE SPIN ASYMMETRY IN THE DRELL-YAN PROCESS 

The final result for the single spin asymmetry in the DY cross section arising from a soft gluon pole contribution 
in a particular twist-three hadronic matrix element |@-Pl as presented in Ref. ||l| is of the form (Eq. (15) in [^1): 
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The unpolarized quark distribution q{x) and the twist-three soft gluon pole function T{x, x) depend on the 
lightcone momentum fraction x of the quark inside the transversely polarized hadron only. The function T{x, x) 
is up to a minus sign identical to the function T{x,St) introduced in Ref. 0. We have chosen to denote the 
mass of the lepton pair by Q, instead of M as was done in ||l|, in order to avoid confusion with a hadron mass. 
The angles are defined in the dilepton center of mass frame, cf. Ref. Q|. 

In Ref. p] a similar result was derived from another perspective, namely that of so-called time-reversal odd 
distribution functions. The asymmetry (Eq. (71) in |g]) was found to beQ (just as Eq. (|l|) given in the dilepton 
center of mass frame) : 
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The first term in the asymmetry (proportional to /t) equals the term proportional to T{x,x) in Eq. (|l|). To be 
more explicit, 
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^We have corrected the proportionality constant, i.e. we have replaced 4 by —2. 



where e;^^" = ^STom+n- _ e^"^'^ STpn+pU^fj and the correlation function <I>^ in the A^ — gauge is given by 
^%,^{x,y) ^ I ^^e^^^e^^^y-^\P,S\-^^{0)F+^ivn.)M^n.)\P,S). (5) 

The hadron momentum is up to a mass term proportional to n+, which is one of two lightlike vectors n+ and 
ri_, chosen such that n^ ■ ri- — 1. We will often refer to the ± components of a vector p, which are defined as 
p = p ■ T7,zp. To be more specific, we have decomposed the hadron momentum and spin vectors P and S as 

5^-^<-^-^- + ^^- (7) 

xMv2 Qv2 

The second term in Eq. d3) is another soft gluon pole contribution to the same single spin asymmetry in 
the DY process. It is not proportional to T{x,x), but to a chiral-odd projection of $^ at the point x = y. 
Apart from this chiral-odd term, the difference between Ai and A2 is in the additional term in Ai proportional 
to xdT{x,x)/dx. To see the discrepancy more clearly we will look at Eq. (^ in case of one flavor, a purely 
transversely polarized hadron, i.e., \Sirp\ — 1, and disregard the chiral-odd term. The definitions of the azimuthal 
angles and (f>Si differ by 7r/2 and the unpolarized quark distribution q{x) = fi{x). Hence, one has for the 
asymmetry 
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which indeed differs from Ai given in Eq. (H) by the term proportional to xdT{x, x)/dx. Such a derivative term 
was indeed found in other processes, like prompt photon production |@| and pion production in proton-proton 
scattering Q, where it arises from collinear expansions of the hard scattering part. As we will show in the 
next section, the derivative term obtained in Ref. [Q for the particular case of the DY process at the tree level, 
arises in a different way, namely from the unnecessary requirement of n_-independence of the Fierz decomposed 
hadron tensor. Of course, such a derivative term can arise from collinear expansions of the hard scattering parts 
beyond tree level. 

Before we discuss this problem, we will first carefully consider the color gauge invariance of the description of 
the cross section, which is an important issue in itself. We will first review the case of deep inelastic scattering 
(DIS), after which we return to the analysis of the DY process. 

III. DEEP INELASTIC SCATTERING 

A. Color gauge invariance 

We will consider the description of the deep inelastic scattering cross section in terms of color gauge invariant 
functions involving nonlocal operators. The hadron tensor PF'^" can be written as pO|,pT|] 



W^'^" = [ dxTT[E^'''{x)<^ix)] + fdxidx2TT [^^''(xi,X2V"^$^(a:i,X2)] +••• (9) 

Here, E'^'^ , Ej^'^ and $, ^^ are hard and soft scattering parts, respectively. The hard part E'^'^ is the discontinuity 
of the 7*-quark forward scattering amplitude and Ej^"^ has an additional gluon connected to it. They are IPI 
in the external legs. The color gauge invariant soft parts are defined as 
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where 
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The projector uj" /^ — g" a — n'^n^fj implies that the index of the covariant derivative is either transverse or in 
the minus-direction (the latter contributes at order 1/Q^ and will be neglected here). The above result holds 
for DIS up to order \/Q and is relevant for polarized scattering. It was arrived at from a collinear expansion of 
the starting expression pO|,pT[| 
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without the projector to and with color gauge variant Green's functions 
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The collinear direction n^ is defined by the hadron momentum (up to a hadron mass term) and the parton 
momenta are decomposed as 
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From power counting one sees that the plus-component of <i>^ can start to contribute at order 1. The transverse 
and minus components are order l/Q and 1/Q^ respectively. In a general gauge one writes A" = n" (n_ • A) + 
uj"/3 A^ , such that the term proportional to n" will yield 
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due to the Ward identity |I( 
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Combined with matrix elements containing arbitrary numbers of A'^ gluons (all contributing at the same order, 
which in this case is order 1), this will exponentiate into the path-ordered exponential in $(2:) |jl3-n3|, with a 
straight path along the n_-direction, Eq. (p^). Hence, color gauge invariance of the factorizcd hadron tensor is 
manifest and for instance, a lightcone gauge A+ — can be chosen without problems, reducing the path-ordered 
exponential to unity. 

We note that one starts out with the correlation function $J|, without restriction on a, and then derives an 
expression involving <I>^, where only the transverse components are contributing up to order 1/Q^ corrections. 
The parametrization of <I>^(a;i,a;2) can then be chosen to reflect the fact that only the transverse part is 
contributing. This is the course taken in Ref. [p6[, which uses the parametrization 
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where et^'^ = e'^^^'^n+a'n-js. Obviously, this parametrization depends on n_, nevertheless, one will obtain an n_- 
independent cross section if n_-independence (or equivalently Lorentz invariance) is imposed on the functions 
in the starting expression Eq. (13). We will address this issue of n_-independence and its consequences in the 
next section. 



B. Lorentz invariance and n -independence 

We have shown that choosing a hghtcone gauge n^ ■ A ~ is not required, hence, the vector n_ is not 
(necessarily) associated with the choice of gauge. The vector n_ reflects the choice of basis in which the hard 
scattering process is expressed. The requirement of n_-independence of the final result is then simply the 
requirement of Lorentz invariance. Since the Lorentz covariant tensor W'^'^ is written in terms of traces of 
hard and soft parts after the choice of this basis, n_-independence of the hard and soft parts separately is 
not a requirement. For a different choice of n_, one would obtain a different expression for the hard and soft 
parts and only the trace of the two parts will be 7i_ independent (which wc have explicitly confirmed for the 
calculation of the DY asymmetry in Ref. H)- The reason that the choice of n_ has to be made is because in 
DIS there are only two external momenta, namely P and q, hence not all directions that arise in the calculation 
can be expressed solely in terms of external momenta. In the DY process this can be done using Pi , P2 , <? and 
(Pi X P2) X q, hence one need not resort to an arbitrary choice of n_ in principle. Therefore, in the case of the 
DY process the issue of n_-independence can be avoided altogether. 

The expression for the hadron tensor in Ref. Q| (taken from Ref. ||ll|) differs from Eq. (||) (after insertion 
of the parametrizations of $ and $f)), due to the additional requirement of ri_ -independence of the factorized 
hadron tensor expression Eq. (^. To be explicit, the expression for the hadron tensor in Ref. ||l| (their Eq. (1)) 
is given by (we have replaced the index /i by a and the hadron momentum p by P to avoid confusion) 

+ J dxi dx2 B''{xi,X2)^Tt [£;„(xi, X2)/'7'] ^" + / d^^ ^^^2 B^{xi,X2)^Tt [E^ixi,X2hp] e""^^, (20) 

with eP"^^ = eP^i^'Sf.Py. The functions E and E^ are the hard scattering parts Ei"" and E^" of Eq. (|). The 
Dirac projections of the soft parts <l>(a;) and $^(a;,y) are denoted by the functions Ct,B^ and B^ and hence 
constitute the parametrization of $2). The function B^ is related to the function Go of Eq. (19) as (cf. Eq. (4) 

ofi) 

B^(xi,X2) = -iTr[$Da(a;2,xi)^_]e"^"+"-/P+ = 2M52Gc(x2,xi). (21) 

The function B^ is related to the function Gd and the chiral-odd functions Hd and Ed are not considered in 
Eq. (§). 

The last term in Eq. (GG) is the one that is relevant for our discussion, since the soft gluon poles reside in the 
function B^ . The relation to T{x,x) is given by 

T(x, x) = lim Tix, y) = lim [ttB^ {x, y) {x - y)] /g, (22) 

which shows that if T{x,sc) is nonzero, the function B^ {x, y) has to have a pole at the point x = y. 
The last term in Eq. (M) is accompanied by the term eP°'SP ^ ^paSn+ ^ jsjgxt we note the identity 



The second and third term on the r.h.s. restrict the index a to be transverse, unlike the first term, which 
implies that a can be +. The first term potentially produces a contribution in the hadron tensor Eq. ( pO| ) that 
is proportional to B^ and hence Gd, whereas the projector lo in the hadron tensor Eq. (g) prohibits such a 
contribution. 

The expression for the hadron tensor Eq. (^0|) can be written in the form of Eq. (|9|) without w projector, in 
which case <&£) is given by 
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which is equivalent to imposing rt_ -independence on the parametrization of $2) (without the w projector), which 
is clearly different from Eq. (|l^) (disregarding the chiral odd pieces) . As said before this will potentially imply 



the appearance of additional contributions. In DIS such additional contributions happen not to ariscQ but in 
the DY process it will be the case. So the two expressions for the hadron tensor (Eq. (Q) and Eq. (po|)) give the 
same results in DIS, but they do yield different results in the DY process calculation. In fact, the first term in 
Eq. ( |23| ) is the origin of the x dT{x, x)/dx term in the asymmetry expression Ai Eq. (|l|). 

In conclusion, the requirement of n_-independence of the hadron tensor Eq. (Bh will lead to a different answer 
for the Drell-Yan cross section in the case soft gluon poles are present. Without imposing this unnecessary 
requirement the derivative term will not be present in the asymmetry expression. Let us emphasize that this 
does not apply to the derivative terms obtained from coUinear expansions of the hard parts, such as arising 
in prompt photon production 0, pion production in proton-proton scattering n and (presumably) the DY 
process beyond tree level. 

IV. COLOR GAUGE INVARIANCE IN THE DRELL-YAN PROCESS 



In section III A on DIS we discussed how the collinear polarization part of any number of gluons (the A+ 
gluons) exponentiate using a Ward identity rendering the nonlocal operators color gauge invariant. In analogy 
to the starting point, Eq. (|l3|), for DIS, for the Drell-Yan process the starting expression for the hadron tensor 
integrated over the transverse momentum qt of the lepton pair and including order l/Q contributions, is given 

by 
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dyTr ('$^(y,a;)7''$(x)7„-^ ^— ^7^) + [ dyTr fe(x,y)7^-^ "" ~ ^ . 7c.Hx)y) 



QV2x-y + ie 'J J ^ y av :«y/ Qy/2 x - y - ie 
dyTr ( ^x)r^"Ay,^h''^ -''Z^ . 7a') - f dyTr fj>(x)7a-^ _ "" 7 ^ . rKi^.vh"' 



QV2x-y + ie J J V QV2 x - y 

In this expression the terms with fi± arise from the fermion propagators in the hard part neglecting contributions 
that will appear suppressed by l/Q^ (for details sec Ref. ||]|j). In case one assumes that the functions ^^ix, y) 
(and $^(a;,y)) to be regular in the point x — y, where the gluon has zero momentum, one can replace terms 
like 

fi+ x-y fi+ ^26) 



V2Qx~y + ie ^Q 



In se ctions |IV A| and |IVB| we will assume that this is indeed the case and we will return to this issue in section 
[VC , when we investigate the consequences of poles in <I>^(2/,a;) at the point x = y, the so-called soft gluon 
poles. 

Below we explain how Eq. ( pq ) can be expressed in terms of manifest color gauge invariant matrix elements 
by summing collinear gluons. One conclusion will be that, also with the effects of transverse momenta included, 
color gauge invariance implies that the index on <I>^ is in fact transverse. 

The issue of color gauge invariance in processes with two soft parts, like the DY process, is a theoretical topic 
not yet fully addressed in the literature. One has to show that a process with two soft parts factorizes into color 
gauge invariant objects. We note that if the cross section factorizes into color gauge invariant functions, this 
will allow one to use different gauges for different correlation functions in the same diagram. 

For DIS, a process with one soft part. Ward identities applied to correlation functions with arbitrary numbers 
of A+-gluons, yield the desired path-ordered exponentials or link operators, that render the correlation functions 



^This is in fact due to time reversal invariance. The derivation of Eq. (EOl) is given in Ref. pil and requires n_ 
independence of the Fierz decomposed hadron tensor. The n_-dependent terms are required to cancel each other, 
but the resulting equation (Eq. (43b) in Ref. uW} involving the function B^ turns out to vanish due to time reversal 
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Note that two misprinted signs for the ie pole prescriptions as given in Ref. pi have been corrected, the results there 
were produced with the signs given here. 



color gauge invariant (see the previous section). In [|12|,y_7 15 1 this issue was considered for the DY process, but 
the transverse momentum of the quarks is not included in the most general way. Unlike [[l2yi7|Jl5[| we will focus 
on the case that the transverse momentum of the lepton pair is not integrated over and is small compared to 
the hard scale. In Ref. gjl^l transverse momentum is taken into account, but only for the leading twist. We 
will first consider tree level in the hard scattering part, but all orders in the coupling constant appearing in the 
link operators. Next we will discuss the extension of these results to all orders in as in the hard scattering part 
by means of Ward identities. Our analysis will show explicitly how coUinear Ward identities can be applied to 
the case of the DY process with transverse momentum, including next-to- leading twist, in order to show that 
coUinear gluons exponentiate. This will also be an important ingredient in a future proof of exponentiation of 
soft gluons into Sudakov factors. 



A. Color gauge invariant correlation functions at tree level 



We will now show that the DY process including l/Q contributions, actually consists of color gauge invariant 
correlation functions containing link operators with straight paths along a lightlike direction extending to 
lightcone infinity. Such a link operator in $(p) is for each colored field of the form 
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(27) 
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We will extensively study the tree level situation and note that the coupling constants appearing in the link 
operators are part of the color gauge invariant definition of the correlation functions, hence, count as tree level 
objects and not as part of the perturbative corrections to the hard scattering. We will consider the case of the 
DY process, but the case of two-hadron production in electron-positron annihilation or semi-inclusive DIS will 
be completely analogous, although the direction in which the links are running is different. For instance, in the 
case of leptoproduction one finds for each field in the correlation function $(p) the link operators 



£^[oo,a;] = Pexp I ig / dy A'^{y) 



(28) 
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Since the cross section as a whole is gauge invariant, we will choose one lightcone gauge A~ = and show that 
summation of A'^ gluons leads to appropriate link operators in one set of correlation functions (the correlation 
functions <&, $f)). This results in the color gauge invariance of these correlation functions and thereby we will 
have demonstrated that choosing a lightcone gauge in correlation functions is not necessary (any or no gauge 
choice will do). 

In the A^ = gauge one has to include y4+ gluons in all possible ways. We observe that y4+ gluons emanating 
from $ will give rise to \/Q^ suppression, so we only need to consider the cases where the A+ gluons emanate 
from $. We will show that matrix elements with multiple A+-gluon fields in <I>^, $^^, . ■ • and $^^, $^^^: • • ■ 
will combine into link operators in $ and $^, with paths as given above (running along the minus direction 
extending to y^ — — oo). Note that the correlation function <I>^ appears in the end, since the explicit A" in <I>^ 
is not gauge invariant. We also note that in the A+ = gauge one has (due to iD°' = id"' + gA°') 



$s(p) ^ l^v n{pi,p) ^ Jd'pi $sbi,p) -p" hp), 



(29) 



but for color gauge invariant correlation functions this relation will be modified, since id" also acts on link 
operators. 

We will first consider the case where all gluons connect to the left side of the diagram. In Fig. n^ we depict one 
such diagram that is relevant for the tree level calculation including l/Q corrections [g. The diagrams with first 
of all zero or one A'^ gluon(s) that we need to consider are (schematically) displayed in Fig. 0. We have only 
drawn the relevant left part of the diagrams and left out the blobs denoting the correlation functions, so gluons 
going to the bottom (top) emanate from $ ($). For instance. Fig. ||b represents Fig. |l|. The index a in Fig. ||b 
can be either transverse (T) or -I- and it turns out that a and /3 in Fig. gc will be T and -I- respectively, otherwise 
the diagram will give a suppressed contribution; on the other diagrams we have indicated the appropriate gluon 
polarization (coUinear or transverse). 

The diagrams (h) and (i) are the only diagrams with a triple gluon vertex, since others will either be of 
higher order in as or cannot be separated from the soft part, since not all momenta are participating in the 
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FIG. 1. A diagram contributing to the Drell-Yan process at order 1/Q. 



hard scattering and soft momentum loops will appear. We will first ignore the issue of color matrices in the 
diagrams without a triple gluon vertex and return to this point after we have considered the diagrams (a)-(g), 
i.e., after Eq. @). 

To discuss the expressions for the diagrams as given in the figure, we will introduce the shorthand notation: 
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Also, we need the precise definitions of $|^(pi,p) and <i>^^ (p2 , Pi , p) {ct,(3 not restricted to be transverse) and 
the antiquark correlation functions $(fc), <i>^(A:i, /c): 
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The hadron momenta Pi, P2 (chosen such that Pi ■ P2 ^ ©(Q^)) and photon momentum q are decomposed as 
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for Qy = — g^ ^ Q . Furthermore, we decompose the parton momenta as (i = 1, 2, . . .) 
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FIG. 2. A subset of diagrams to be considered upon inclusion of zero or one A'^ gluon(s). 






We find for diagrams (a)-(d): 
(a) = f$(p)7''$(fc)7'^ 



(&)=||<i'Sb)7^'J>(fc)7aT^7'^ 



Ai^ibi'P) 



x/2 



Qil-^i+ie, 



■7^$(fc)7'- 



dV *i(pi'P)^''*(fc) Q(r^f:p-) {^^^^ + ^^- ^1^) - ^- h'' 



(c) 



V2 



^^P. Ar a>-(,.p„p) ^^^-^^^ 7'^<P(fc)7..^7^ 



J4 ,„ 7S-/3 , 



(d) = / dV$+(pi,p)7^ / d^fci<i>^(fci,fc)7;3T + $(fc)^: 



V2 



0(1 - 6 + ie) V2Q 



(41) 



(42) 



(43) 
(44) 
(45) 



In the last equation we have used the fact that in the A = gauge no link operator will be present in $^. 
We can then apply the equations of motion to j d'^ki $2?(^ii k)^pT^ which results (in the A^ = gauge) in a 



cancellation of the term with the P_ = fi-fL+ji projector in the expression for diagram (b). Anticipating the 
specific link operators that will result in the end, we observe that 



$[-oo,a;](p) 



■)- 
d'^x 



72 



$(p) + y'dV$i(Pi,p)Q^Y:^ 



+ ie) 



(46) 



(2^)4 



V2 



•I'aW 



y d^P2 d^pi $^+ (p2 , Pi , p) Q^^_^^ 



(47) 



where the definition of $^(p) is given in Eq. (P9|). The smxi of the four diagrams (a)-(d) therefore results in (to 
order g in the link operator) 



(a) + (b) + (c) + (d) = ||$[-^,a:](p)7'^<i>(fc)7l| + || $Shcx3, x](p)7''$(fc)7aT^7'' 



dV$i(Pi,p)7^$(fc) 



V2 yi+ 

Q(l - Ci + le) V2Q 



HT 



)Y 



(48) 



By taking into account the derivative of the link operator, the last two terms nicely combine into the following 
expression containing exclusively the functions <&[— oo,x](p) and ^'^[—oo,x\{p) (the latter is obtained from Eq. 
(|^) by replacing gA" by iZ?"): 

(a) + (&) + (c) + {d) = |$[-cx3,a;](p)7^$(fc)7'' 



{$^ [-00, x] (p) - p" $[-cx), a;] (p)} 7^$(A;)7aT 



The diagrams (e)-(g) also turn out to combine, yielding the (expected) expression: 



^/2Q 



(e) + (/) + {g) = -[$[-oo,x](p)7^ {7o<(fc)7o + fc'' <&(&)} t'^^T/st 



(49) 



(50) 



74),.7f:^ , 



where we used that J d'^ki^^jj{ki,k) = 7o<&_d (fc)7o, when <i>£,(fc) = J d^ki^jj{k,ki), cf. Rcf. ||T^. To arrive at 
Eq. (pG) we also made use of the fact that to this order $^ (with a transverse) will be proportional to fi- and 
that $^ will be proportional to ^+. 

Until now we have ignored the issue of color matrices. It is easy to see that in diagrams (d), (f) and (g) the 
color matrices are in the wrong order for them to be absorbed into the correlation functions along with the link 
operators. The diagrams with the triple gluon vertex should be included to cancel the contributions arising 
from commuting the two color matrices. In this way diagram (h) will cancel the term arising from diagram (d) 
and diagram (i) that of diagrams (f) and (g). Hence, by including diagram (h) and (i) one can include the color 
matrices in the link operators. 

This concludes the case of including zero or one A^ gluon(s) (emanating from <&, connecting to the left). 
Inclusion of two A"*" gluons is done in a similar fashion, where we note that from diagram (c) with a = /3 = + 
one can see that a product of two theta-functions in the minus component appears, which gives rise to the path- 
ordering. Knowing this, it is a matter of iteratively treating each order like the first order we have discussed 
just now, by considering {$g[— oo,a;](p) — p" $[— oo,x](p)} as an effective $^(p). 

When the A^ gluons emanating from <I> connect to the right side of the diagrams, then one arrives at the 
link operator 



£+ [0, -(^] = P exp ( ig J dy- A+ {y) 



(51) 



y+=0+,yj,=Oj 
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Inclusion of arbitrary numbers of A"*" gluons emanating from <I> connecting to either side of the diagram, results 
in the appearance of the following color gauge invariant correlation functions exclusively: 



$Gi(p) 



{2n)' 



e'P-^ (Pi, 5i|^(0) £+[0, -(»] £+[-00, x] 7A(a;)|Pi, 5i), 



(52) 



'^ ^ ^ y „wy „wi-(^-v) 



(Pi, ^ilV(O) £+[0,-00] 



(27r)4 (2^)4 
X C+[-^, y] iD'^iy) C+[y, -00] £+hoo,x] V(a;)|Pi, ^1). 



(53) 

Hence, Eqs. (p2[) and ( p3| ) are the Drell-Yan process extentions of Eqs. (Eoh and (0) to the case the correlation 
functions not only depend on the lightcone momentum fractions. 

So we conclude that at the tree level the complete set of diagrams with arbitrary numbers of A+ gluons 
emanating from $ will result in the appearance of color gauge invariant correlation functions (containing link 
operators and covariant derivatives), with path-ordered exponentials having paths along lightlike directions 
extending to lightcone infinity. Not choosing the A~ = gauge cannot affect $gi or <&^qi to this order (1/Q) 
and hence, these are the correlation functions that will appear in any gauge. Choosing the A~^ — gauge 
instead, will give rise to similar results for $ and $£,. The color gauge invariant matrix elements <I>gi(^) and 
^jjQi(ki,k) contain similar links extending to lightcone infinity, but now involving A~ gluon fields: 

¥Gi(fc)= /44je-'''"(^2-^2|^(a:)£^[x,-oo]£-[-oo,0]^(0)|P2,52), (54) 



(2^)4 



$ 



DG1\ 



(fci,fc)^ /_^_^e-H.-.)e-^-.( 



(2.)M27r)4^ ^ ,P2,52|^(x)/:-[x,-oo] 

X £-[-00, y] ii?"(2;) £-[2/, -(^] /:-[-(^, 0]^(0)|P2, ^2), 



where 



C [x, —00] = V exp ig j 

\ J —00 



dy+A (y) 



(55) 



(56) 



/ y =x- ,yT=XT 

One thus finds the following color gauge invariant expression for the hadron tensor not integrated over the 
transverse momentum of the lepton pair, including 1/Q contributions: 



W''''{Pl,P2,q) = ^ J d^kd'pS\p + k-q)lTT[^Gl{p)r^Gl{kh''] 



Tr 



Tr 



Tr 



Tr 



Al <^DGl{Pl,P)~P°' *GI(P) } 7^*Gl(fc)7aT 



V2Q 



^Gi{p)Y' [Jd^ki <^DGiiki,k) + kP <l>Gl(fc)| I'^IPT 



rfV$SGl(P,Pl) -P"<&Gl(p) [>7^-^7aT$Gl(fc)7'^ 



^GMltiT^r \ I d^k, <Gi(fc, fci) + fc'^ $Gl(fc) I' Y 



(57) 



This expression still contains some order 1/Q^ contributions, which we ignore to keep the notation as simple as 
possible. To the order we consider one can for instance replace 

(S'*(j9 + k-q)^ 5{q+ - p+) 5{q- - fc") 5^{pj, + fcj. - qj,). (58) 

To obtain the hadron tensor expression Eq. (|5^) we have already used these delta functions, for instance in the 
parametrizations of the momenta g,p, k. 

We would like to emphasize that we have not obtained color gauge invariance at the cost of introducing a 
path dependence. The expansion in orders of l/Q does not leave any freedom of choosing the paths for the 
leading and next-to-leading twist terms. 
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B. Color gauge invariant correlation functions to all orders in as 

In order to extend the above tree level analysis to all orders in as in the hard scattering part, we will 
first discuss two limiting cases. The first case we consider is the cross section integrated over the transverse 
momentum qrp of the lepton pair. The second case is the cross section at measured q^ (with q^ = Q\ ^ Q^), 
but restricted to the leading twist. 

In case one integrates the cross section, and hence the hadron tensor, over the transverse momentum Qj, of 
the lepton pair, then the situation simplifies considerably. After a coUinear expansion of the hard scattering 
parts, the transverse momenta of the partons can be integrated over each separately. The hadron tensor will 
be expressed in terms of coUinear partons only, i.e., they are collinear to the momentum direction of the parent 
hadron (neglecting target mass corrections). More specifically, the hadron tensor contains the above gauge 
invariant correlation functions $, $^ ($, $£,) integrated over the — (+) and transverse momentum components. 
They depend only on the lightcone momentum fractions, such as the ones that appear in the DIS cross section, 
cf . Eq. (|l^ and ( |ll| ) . In that case one can apply Ward identities to gluon fields with polarizations proportional 
to their momenta, to prove that the above results are also correct to all orders in as in the hard scattering part. 
This has been done for the leading twist contribution in for instance Refs. [n3,E0|, Ref. [O (cf. its Fig. 4.4 and 
appendix) and Re f. p4| (cf. its Figs. 19 and 25 and Eqs. (131) and (133)). It is completely analogous to the 
application of Eq. (|1^) in DIS. Also, it straightforwardly applies to the 1/Q corrections for which a factorization 
theorem has been discussed in Ref. [£1[. The same Ward identities can be applied, since the additional At 
gluon has physical polarization and is on-shell and therefore, does not affect how the collinear gluons appear in 
the Ward identities. The latter just produce link operators on both side of the At field and since in this case 
the links are independent of transverse coordinates, the replacement qAt — > iDt — iOt immediately gives the 
result for J (Pqrp)/V^''{Pi,P2,q), using VV^"^ from Eq. ( p7\ j with hard scattering parts H containing all orders 
in 0;^: 

J d^QT W'^-'ix, X, <It) = \\ Tr [$Gi(a;)$Gi(S)i?^"(x, x)] 

dy {•^'hG,{y,x)-5{y~x) [ d^pS{p+-xP+)p'^'PGi{p)}^Gi{x)HZ{y,x,x) 



-Tr 



-Tr 



$Gi(x) J dy U^'^G.iy^x) + 6{y - x) J d^k6{k- - xP^)k^<PGi{k)\ HZ{x,y,x) 



(59) 



This is just a schematic way of writing, since the four spinor indices on the hard scattering part iJ'"^ and Hl^'^ 
connects to both $ and $. At tree level this equation is Eq. (p5|) re-expressed in terms of color gauge invariant 
functions, Eqs. (|lO|), (|Il|) and the obvious extensions $Gi(a;) and $^ gi(^' v)- This holds in case soft gluon poles 
are assumed to be absent. Our previous tree level analysis is in this case an explicit demonstration of a result 
valid to all orders in a^, since soft gluons can be shown to cancel pi|| . 

To apply such Ward identities to the more extended case we considered here, namely to the cross section 
differential in the transverse momentum q^ of the lepton pair (Q^ ^ Q^), is more complicated. The momenta of 
the partons are in general not collinear to the parent hadron, since they also possess some transverse momentum. 
Hence, the A+ gluons, i.e., the gluons that are collinear to the parent hadron do not have polarization collinear 
to their own momenta. The latter type of gluons are called gluons with longitudinal polarization mm and are 
the gluons to which Ward identities can be applied directly. For the leading twist contribution the application of 
such Ward identities still works for ^4+ gluons, because deviations of A~^ gluons from longitudinal polarization 
(such deviations are proportional to At) will be suppressed by inverse powers of the hard scale. We then find 
for the leading twist part of Eq. (|5^) with a hard part H containing all orders in ag : 

W-^ix, X, Qt)^\ I d^PT d^kj, S^{Pt + ^t " Qt) Tr [^gi{x, Pt) *gi(x, kr) il^"(x, x,Pt, k^, Qt)] ■ (60) 

However, this is not the complete result, since if one goes beyond the tree level, then apart from having a more 
complicated hard scattering part, soft gluons have to be resummed into Sudakov form factors B, resulting in 
a replacement 

S^ip^ + k^-q^)^J^ e-^HPr+'^r-lr) ^-^if) ^ (61) 
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where e ^'■■^^ is the Sudakov form factor. This has been shown in Ref. [g 18 1 for the leading twist 



For the next-to- leading twist a factorization proof for the cross section differential in the transverse momentum of 
the lepton pair {Q"^ ^ Q^) has not been given. But we will discuss how Ward identities applied to longitudinally 
polarized gluons can be used to show that coUinear gluons factorize, i.e., they exponentiate into link operators, 
also at next-to-leading twist {0{l/Q)). 

From our tree level result we can deduce how the deviations of longitudinally polarized gluons from coUinearity 
will affect the application of the Ward identities. The part of a longitudinally polarized gluon field A^ ~ p^ 
that deviates from the n+ direction (i.e., from its A'^ part) will be proportional to its transverse momentum 
(^ipT ^ Pt)- ^^'^ latter part will generate derivatives of the A'^ parts of the other additional longitudinally 
polarized gluons, which to this order can be taken to be coUinear. This yields the derivative of a link operator 
that is exactly needed to compensate for the effect of the replacement gA -^ iD — id inside ^'\ (which now 
contains link operators depending on transverse coordinates), in order to arrive at the replacement <&^(p) -^ 
^%{p) — p"^(p) for correlation functions containing link operators. At tree level this is the step from Eq. (H^) 
to Eq. dl). 

To make this a bit more explicit we write down the longitudinally polarized gluon field in Fig. y 



^rp = n+(^--^) + (P?-P?T) 



V2 



Qil-Ci+ie) 



{n.-A) ex (p"-p?). 



(62) 



Insertion of this field into the most general hard part (Fig. n^ contains only the tree level hard part) and applying 
the Ward identity will result in those diagrams where the longitudinally polarized gluon attaches directly to the 
"external" legs of the hard part (those legs that connect to the soft parts), cf. Ref. pj. The v4+ part of A? 
will yield a term in a link operator in <& (represented by a double line in Ref. ||lj], called an eikonal line) and 
the second term on the r.h.s. of Eq. (p2h will give rise to the derivative of a link operator in ^'^. At tree level 
these are the first and third term of Eq. (|4^) , respectively. 

It is easy to verify that this will work to all orders in a^ in the hard part, since one has to take this deviation 
into account for only one longitudinally polarized gluon. Hence, our conclusion is that transverse momentum 
dependence in correlation functions does not spoil the summation of coUinear A"*" gluons into link operators at 
next-to-leading twist (0{1/Q)). This is not in conflict with the known failure of factorization at even higher 
twist |22). 

Hence, Eq. dS^) can then be written with a hard part H containing all orders in as: 



W''''iPi,P2,q) = ^ J d^kd^pd\p + k-q)lTT[<^Glip)'^Gl{k)H'^''ip,k,q)] 



Tr 



Oi {nGiiPi^P)-SiPi-p)p"^Gi{p)} >^Gi{k)HZ{p,k,q) 



Tr 



*Gi(p) 



d^fci 



{<Gi(fci, k) + 6{ki - k) k" $Gi(fc)} H^^ip, k, q) 



(63) 



Again soft gluons will have to be included and they will modify the delta function. 

Application of Ward identities are an essential part of the proof of factorization. They are used to show that 
coUinear gluons exponentiate, like discussed above, but also that soft gluons factorize from the hard part and 
from the correlation functions. In the grp-integTaAed case the soft gluons must be shown to cancel and in the 
case of measured g^ they must be shown to exponentiate into Sudakov factors. For the next-to-leading twist 
the latter case has not been studied yet, but the application of coUinear Ward identities as discussed above 
will be an essential ingredient. The fact that Ward identities can be used to show that coUinear gluons also 
exponentiate at next-to-leading twist in the case of measured qrp (with Q^ -^ Q^), suggests that coUinear Ward 
identities can also be applied to the soft gluon exponentiation and a factorization theorem is expected to hold. 
It could of course be the case that the replacement Eq. (|6l|) is different for these subleading terms than for the 
leading term. This has to be investigated further. 

Finally, we have to remark that additional divergences are introduced by the paths extending to infinity if 
one considers as corrections connecting to the paths; these have to be renormalized. A discussion on this topic 
and how to deal with it is given in 0] and applies unchanged to the next-to-leading twist results discussed here. 
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C. Color gauge invariant soft gluon poles 

Until now we have considered the case where we assume all fields at infinity to vanish inside the hadronic 
matrix elements. In the case of soft gluon poles this assumption is relaxed in a specific way. To include such 
contributions one has to make the following replacement (cf. Eq. (^5|)) in the first term in Eq. (|4^) for diagram 
(b) and in Eq. (Q for diagram (c): 

^+ , ^+ 1-^1 (64) 



\/2Q \/2gi-6+ie' 

A soft gluon pole means that the correlation fmiction has a pole at ^i = 1, i.e., when the gluon has vanishing 
momentum. This prohibits applying the inverse of the above replacement . 

If we consider the above replacement the following matrix element (times (5(p^ ^ Pi )) results (in addition to 
the matrix elements in Eq. (pTj)): 



(27r)4 (27r) 



|y e*P-y ^^Pi-i'-c-y) (p^^ 5-^1^(0) £+[0, - 



-00 



x\^{C+[-^,y]gA^{y)C+[y,-^]) C+[-^,x]i^{x)\P^,S,). (65) 

Due to time reversal symmetry and under the assumption of continuity of the correlation functions, one can show 
g that for soft gluon poles one has to have antisymmetric boundary conditions on the field At at y^ = ±c». 
In this case one can replace gAx -^ iDt in the above matrix element. The identity 



_dy 
is used to arrive from Eq. (pq) at 



^-,C+[-^,y]iD'-{y) C+[y,-^] 



£+h(X3,2;]5F+"(y)/:+[y,-(x.] (66) 



a>^Gi(Pi, P) = /(^(^e*^-^e^''-(--^) (Pi, ^ilV'(O) £+[0,-00] 



xC+[-^,y]^-F+'-{y)C+[y,~^]C+[-^,x]i;{x)\P,,S{), (67) 



2 

where a is transverse. The gauge invariant term $^gi(pi,p) will show up in the hadron tensor Eq. (pTh in 
exactly the same way as $^Qj(pi,p). In Eq. ( p7| ) one can just replace 

^iGliPl.P) ^ ^'hGl{Pl,P) + nGliPuP)5{p+-pt). (68) 

We find that upon choosing the A+ — gauge in J d'^pid'^p^Q,j^^{pi,p)5{p^ ^pt) ^ti^^ the standard expression 
for the soft gluon pole matrix element T(x, St) as given in ^,||| results, cf. Eqs. (^) and (H) (note that we have 
now included a factor g/2 in the matrix element ^pQi)- 

From the above it is clear that soft gluon poles require nonzero matrix elements containing At fields at 
y~ = ±cx) in any gauge. On the other hand, we have chosen hadronic matrix elements with A'^ fields at 
y~ = ±00 to vanish. The reason is that otherwise the previously obtained matrix elements will not be color 
gauge invariant anymore and their derivation will not even be well-defined, i.e., will not result in finite answers, 
since 

will be divergent. We simply conclude that soft gluon poles arise from gluon fields with physical (transverse) 
polarizations. 

Hence, the inclusion of soft gluon poles can be done in a color gauge invariant way, even in the case where 
the transverse momentum of the lepton pair is not integrated over. 

As a final point we note that in the color gauge invariant correlation functions ^^ qj and ^p qj the index 
a is purely transverse, reflected in the appearance of 7^ in Eq. (^^ (analogous to the projector u) in DIS). 
Imposing the additional, but unnecessary constraint of n_ -independence at this stage will alter the expression 
for the hadron tensor and may result in different answers. 
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V. CONCLUSIONS 

We have considered the color gauge invariancc of a factorized description of the Drell-Yan process cross 
section. The analysis focused on next-to-leading twist contributions for polarized scattering {0{1/Q)) and on 
the cross section differential in the transverse momentum of the lepton pair {Q\ <C Q^). The hadron tensor 
has been expressed in terms of manifestly color gauge invariant, nonlocal operator matrix elements. The tree 
level case was worked out explicitly and Ward identities, applied to longitudinally polarized gluons, allowed to 
extend the results to all orders in the hard scattering part. This use of Ward identities in case longitudinally 
polarized gluons are not equal to coUinear gluons will also be an important ingredient in a full (and still lacking) 
factorization proof for the case under consideration, to be given in a gauge independent way. The factorization, 
i.e., exponentiation into Sudakov factors, of soft gluons is expected to be analogous to the factorization and 
exponentiation of collinear gluons as derived here. 

The leading twist part of the expressions we have derived, are in agreement with earlier investigations pJlSl. 
Moreover, the whole result (including the next-to-leading twist) reduces to known results [ pl[ upon integration 
over the transverse momentum of the lepton pair. We have also given a color gauge invariant treatment of soft 
gluon poles, showing that in any gauge such poles arise from gluon fields with physical (transverse) polarizations. 
In addition, we have extensively discussed the discrepancy between the results of Ref. ||l| and |^ for a single 
transverse spin asymmetry in the Drell-Yan process. This asymmetry appears at the tree level in case soft 
gluon poles in the twist-three matrix elements are present. We have demonstrated that the requirement of n_- 
independence of the resulting factorized hadron tensor expression as imposed in Refs. ||l|,|l[ is unnecessary from 
the point of Lorentz and gauge invariance and leads to a different answer in the case of the above mentioned 
asymmetry. Without imposing this requirement the derivative term will not be present in the asymmetry 
expression at tree level. We have carefully considered the color gauge invariance of the description of the cross 
section in terms of the matrix elements, to conclude that only physical, i.e., transverse gluon degrees of freedom 
contribute and hence, arrive at the expression A^ (Eq. (||) or equivalently Eq. (^)) for the asymmetry. 
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